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ABSTRACT 


A  numerical  lifting  surface  procedure  is  developed  specifically  for 
flapped  control  surfaces  with  trapezoidal  planforms.  The  procedure  uses 
a  discrete  vortex  approximation  with  spanwise  vortex  lines  located  at 
constant  percentages  of  the  chord.  Use  cf  the  procedure  in  obtaining  an 
optimum  flapped  rudder  design  is  demonstrated.  A  listing  and  user’s 
description  of  the  computer  program  is  included. 
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NOMENCLATURE 

(In  the  order  of  appearance  in  the  text) 


a 

s 

A 

f 

Af 

X 

®T 

CR 

A 

x 

z 

XL(z) 

X1,(z) 

i 

j 

m 

n 

XFR 

XFT 

XSR 

XST 

iv 

Ih 

I 

J 

NF 

NS 

G(x,z) 

a 

6 


(2s)2 

effective  aspect  ratio  — 

semi-span 
rudder  area 
flap  area  ratio  A^/A 
flap  area 

taper  ratio  c^/c^ 
tip  chord 
root  chord 

sweep  angle  of  1/4  chord 

chordwise  coordinate  axis 

spanwise  coordinate  axis 

leading  edge  chordwise  coordinate 

trailing  edge  chordwise  coordinate 

spanwise  panel  index  of  the  lattice 

chordwise  panel  index  of  the  lattice 

spanwise  index  of  control  points 

chordwise  index  of  control  points 

flap  chord  at  the  rudder  root 

flap  chord  at  the  rudder  tip 

skeg  chord  at  the  rudder  root 

skeg  chord  at  the  rudder  tip 

spanwise  precision  number 

chordwise  precision  number 

number  of  chordwise  panels  on  rudder 

number  of  spanwise  panels  on  rudder 

number  of  spanwise  panels  on  flap 

number  of  spanwise  panels  on  skeg 

nondimensional  circulation  distribution 

rudder  angle  of  attack 

flap  deflection  relative  to  skeg 


NOMENCLATURE 

(cont.) 


Y(x,s}  dimej sional  circulation  distribution 

U  velocity  at  infinity 

c^£  node  amplitudes 

k  spanwise  node  index 

&  chordwise  mode  index 

K  number  of  spanwise  modes 

L  number  of  chordwise  modes 

fj^z)  kth  apanwise  mode 

■n, 

p^(s)  £,th  chordwise  mode 

r\j 

z  transformed  spanwise  coordinate 

s  transformed  chordwise  coordinate  on  rudder 

c(z)  rudder  local  chord 

r\j 

t  transformed  chordwise  coordinate  on  flap 

Pjj  integral  of  £th  chordwise  mode 

r(z)  spanwise  circulation  distribution 

5,5  coordinates  of  a  general  point  on  a  vortex 

v^  velocity  induced  by  spanwise  vortex 

vCt)  velocity  induced  by  trailing  vortex 

vm,n,i,j  velocity  induced  by  the  (i,j)th  element  at  the  (m,n)th 
3  control  points 

vm  n  k  1  velocity  induced  at  the  (m,n)th  control  point  by  the 
3  *  *  (k,l)th  mode  of  unit  amplitude 

Vm  n  velocity  at  the  (m,n)th  control  point 

p  mass  density 


(^(z) 

CL<$<Z> 

cla 

CL6 

cDi 

T) 


local  lift  coefficient  per  unit  angle  of  attack 
local  lift  coefficient  per  unit  flap  deflection  angle 
overall  lift  coefficient  per  unit  angle  of  attack 
overall  lift  coefficient  per  unit  flap  deflection  angle 
induced  drag  coefficient 
lifting  surface  efficiency 


% 

xH(z) 

c(z) 


first  moment  of  the  Hth  chordwise  mode  function 

local  distance  of  the  center  of  pressure  as  a  fraction 
of  the  local  chord  from  the  hinge  line 
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NOMENCLATURE 

(cont.) 


c(z) 

*h 


°DV 


% 


local  distance  of  the  center  of  pressure  as  a  fraction 
of  the  local  chord  from  the  leading  edge 

resultant  chord wise  position  of  the  center  of  pressure 
relative  to  the  flap  hinge  line  as  a  fraction  of  the 
mean  chord 

mean  chord  -  (xj^  4-  Xpj  +  xSR  4-  xSj)/2 
viscous  drag  coefficient 

transformed  chordwise  coordinate  evaluated  at  the  posi¬ 
tion  of  the  flap  hinge 
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1.  INTRODUCTION 


This  study  arose  from  the  need  to  develop  a  rational  basis  for  the 
selection  of  optimum  geometric  characteristics  for  rudders  with  relatively 
small  flaps.  In  this  case,  very  minor  changes  in  sweep  or  taper  can  effect 
a  major  change  in  the  spanvise  distribution  of  flap  chord.  This,  in  turn, 
could  be  expected  to  influence  the  spanvise  and  chordwise  distribution  of 
lift  when  the  flap  is  deflected.  To  optimize  rudder  performance  one  would 
like  to  have  a  distribution  of  lift  which  results  in  the  maximum  lift/drag 
ratio,  highest  possible  stall  angle,  and  minimum  control  moment  for  both 
skeg  and  flap.  Consequently,  it  is  necessary  to  have  the  means  for  estimating 
spanwi.se  and  chordwise  distributio.i  of  lift  for  a  given  geometry. 

Ship  rudder  effective  aspect  ratios  typically  fall  in  the  region  where 
neither  high  aspect  ratio  nor  low  aspect  ratio  theories  are  valid.  One  must, 
therefore,  resort  to  numerical  lifting-surface  theory  to  obtain  meaningful 
results . 

It  seemed  most  expedient  for  this  application  to  write  a  specialized 
computer  program  designed  to  accommodate  only  trapezoidal  planforms  of  the 
form  shown  in  Fig.  1.  The  flap  hinge  is  required  to  be  at  right  angles  to 
the  root  section,  and  the  tip  chord  is  required  to  be  parallel  to  the  flow. 

No  restriction  is  placed  on  aspect  ratio,  sweep  or  taper,  provided  that  the 
flap  hinge  emerges  from  the  tip,  rather  than  the  leading  or  trailing  edge. 

A  trapezoidal  planform  makes  a  discrete  vortex  lifting-surface  model 
relatively  simple.  Unlike  the  original  work  of  Faulkner  [1] ,  or  current 
schemes  for  propellers  [2],  the  present  work  employs  spanwtse  vortex  lines 
located  at  constant  percentages  of  the  flap  and  skeg  chord,  rather  than  at 
right  angles  to  the  oncoming  flow.  This  eliminates  the  problem  of  vortex 


Fig,  1  Coordinate  System  and  Notation  for  Plaaform  Geometry 
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elements  running  out  of  the  leading  and  trailing  edges,  which  is  inevitably 
a  source  of  inaccuracy.  This  would  be  particularly  objectionable  in  the 
tip  region  of  a  small  flap. 

As  a  result,  the  spanwise  vortex  lines  are  not  purely  bound  vortices, 
since  they  include  a  component  of  vorticity  parallel  to  the  oncoming  flow. 
However,  as  long  as  the  accompanying  system  of  trailing  vortices  is  arranged 
in  such  a  way  that  continuity  of  vorticity  is  preserved,  this  is  an  equally 
valid,  discrete  representation  of  a  continuous  vortex  system. 

It  was  decided  to  use  a  strictly  linear  lifting -surf ace  theory,  since 
an  examination  of  results  for  two-dimensional  airfoils  with  flaps  [3]  indi¬ 
cates  that  flap  effectiveness  is  reduced  due  to  viscous  effects  before  any 
non-linear  augment  of  lift  becomes  apparent.  Consequently,  the  vortex  system 
is  located  on  a  plane,  even  when  the  flap  is  deflected,  and  the  influence 
of  thickness  on  loading  vanishes.  One  may  then  solve  the  problem  of  a  rudder 
with  zero  flap  deflection  and  unit  angle  of  attack,  and  a  separate  problem 
of  a  rudder  with  zero  angle  of  attack  and  unit  flap  deflection.  The  solution 
to  any  combination  of  angle  of  attack  and  flap  deflection  is  then  simply  a 
linear  combination  of  the  preceding  two  results. 
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2  .  DISCRETE  VORTEX  ARRANGEMENT 

The  computer  program  is  designed  to  accommodate  any  quadrilateral 
flapped  control  surface  with,  a  hinge  axis  normal  to  the  flow,  and  with  a 
tip  parallel  to  the  flow,  as  indicated  in  Fig.  1.  The  planfora  geometry  is 
uniquely  specified  by  the  four  nondimensional  quantities  tabulated  below: 


Symbol  Definition 

(2a)2 

a  Effective  aspect  ratio 

f  Flap  area  ratio  A^/A 

A  Taper  ratio  cT/cR 

A  Sweep  angle  of  1/4-chord 

From  these,  we  may  obtain  the  coordinates  of  the  four  corners  of  the  control 
surface.  The  coordinate  system,  as  shown  in  Fig.  1,  1g  located  with  the 
x-axis  situated  at  the  root  section  and  the  z-axis  coincident  with  the  flap 
hinge  axis.  The  semispan,  *i,  is  taken  to  be  unity,  so  that  all  length 
dimensions  are  nondimensionalized  at  the  outset  with  respect  to  this  qr  intity. 
The  x-coordinates  of  the  four  corner  points  then  become: 


3(1  -  A)  +  4f (1  +  X) 
XFR  *  2a (A  +  1) 


tan  A 


XSR  *  XFR  a(A  +  l) 

=  JO.-  1)  +  —■  tan  A 

2a(A  +  1)  * 


XST  a  (A  +  1) 


4A 


(2.1) 
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If  we  do  not  wish  to  have  the  flap  binge  emerge  from  either  the  leading  or 
trailing  edge,  it  is  necessary  that  the  choice  of  input  quantities  be  such 
and  Xpj,  >  0  and  Xg^  and.  Xg^  <0.  If  these  conditions  are  not  met. 
an  error  message  is  printed. 

The  two  trapezoidal  regions  representing  the  skeg  and  flap  may  now  be 
subdivided  into  a  lattice  of  spanwise  and  trailing  discrete  vortex  lines. 

An  individual  segment  of  a  spanwise  vortex,  together  with  the  two  trailing 
vortices  originating  at  the  enu&  of  the  segment,  forms  a  horseshoe  vortex  of 
constant  strength,  as  illustrated  in  Fig.  2. 

The  complete  lattice  arrangement  is  shown  in  Fig.  3.  The  fineness  of 
the  grid  is  controlled  by  specifying  a  spanwi.se  precision  number,  1^  =  0, 

1,  or  2,  and  a  chordwise  precision  number,  Ig  **  0,  1,  or  2.  A  zero 
precision  number  denotes  the  coarsest  possible  grid  spacing,  which  is  the  one 
illustrated  in  Fig.  3. 

The  semispan  is  divided  into  51^  +  7  chordwise  panels  of  equal  width. 

The  panel  nearest  to  the  tip  is  further  subdivided  into  14  equal  intervals, 
and  this  pattern  is  dupJicated  on  the  image  side  of  the  plauform.  The  total 
number  of  chordwise  strips  over  the  span,  I,  for  each  precision  number  is  as 
follows : 


Precision  Number 


Number  of  Chordwise  Panels 


I 

40 

50 

60 


0 

1 

2 


x«,j  ,Z| 


U —  (j  — ►» 

\ 


control  point 
at  xm,n«zm 


Fig.  2  Illustration  of  a  Typical  Horseshoe  Vortex  Element 
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The  total  chord  is  divided  into  J  spactrf.se  panels  in  accordance  with 
the  specified  horizontal  precision  nuaber  1^, 

J  -  50  +  10  Ifl  .  (2.2) 

In  general*  it  is  not  possible  to  have  the  chordvise  spacing  the  seme  on 
both  the  skeg  and  the  flap.  However)  the  two  spacings  will  be  very  nearly 
equal  if  the  nuaber  of  panels  on  the  flap,  HF,  is  the  integer  closest  to 
the  flap  area  ratio  times  the  total  number  of  intervals  over  the  chord, 

NF  >  f*J  .  (2.3) 

The  number  of  panels  on  the  skeg,  NS,  is  then: 

NS  -  J  -  NF  .  (2.4) 

The  determination  of  NF  and  NS  is  internal  to  the  computer  program,  based 
on  the  specified  input  value  of  f  • 

Velocities  are  computed  at  80  control  points  distributed  over  the  span 
and  chord.  Their  spanwise  placement  is  permanent,  as  indicated  by  small 
circles  in  Fig.  3.  Their  chordwise  placement  is  arbitrary,  and  is  specified 
as  input  data  by  the  user.  Spanwise  vortex  lines  are  placed  in  the  middle 
of  each  panel.  Each  control  point  is  located  midway  between  adjacent  span- 
wise  and  trailing  vortex  elements,  which  is  essential  in  order  for  the  dis¬ 
crete  system  to  converge  to  the  Canchy-Principal  Value  of  the  corresponding 
continuous  singular  integral.  Thus  the  control  points  are  placed  on  the 
panel  boundaries. 

To  avoid  errors  due  to  edge  effects,  control  points  should  be  located 
at  least  2  spaces  away  from  the  leading  and  trailing  edges,  and  4  spaces 
away  from  the  tip.  The  choice  of  8  spanwise  and  10  chordwise  stations 
follows  from  our  experience  with  similar  computing  schemes  for  propellers  [2] 
Since  the  flap,  when  deflected,  acts  to  some  extent  as  an  independent  lifting 


surface,  it.  is  important  to  include  as  many  control  points  as  possible  over 
the  flap  chord,  while  avoiding  the  immediate  vicinity  of  the  hinge  line. 

This  requires  an  unusually  large  number  of  panels  over  the  chord,  and  a 
flexible  system  of  specifying  control  point  locations. 

Tills  lattice  arrangement,  which  concentrates  at  least  half  of  the 
elements  over  the  outer  seventh  of  the  span,  is  the  result  of  considerable 
numerical  experimentation.  Lifting  surfaces  with  trapezoidal  planforms  may, 
und»r  certain  circumferences,  have  a  spanwise  load  distribution  which  falls 
to  zero  very  abruptly  at  the  tip.  To  obtain  accurate  results  it  is  necesrjary 
to  locate  a  sufficient  number  of  control  points  very  close  to  :he  tip. 
However,  to  avoid  edge  effects,  as  we  have  noted  before,  it  is  essential  to 
have  several  grid  spaces  between  the  outboard  control  point  and  the  tip.  One 
must,  therefore,  have  an  extremely  fine  grid  in  this  region.  Since  a 
relatively  coarse  spacing  provides  ample  accuracy  in  the  midspan  region  of 
the  lifting  surface,  it  would  be  an  extreme  waste  of  computer  time  to  extend 
the  fine  grii  over  the  entire  planform. 


3.  CONTINUOUS  AND  DISCRETE  CIRCULATION  DISTRIBUTION 
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One  could,  in  principle,  3olve  for  the  circulation  of  each  discrete 
spanwise  vortex  element  necessary  to  satisfy  the  flow  tangency  boundary 
condition  at  each  control  point.  However,  this  would  require  that  the 
number  of  control  points  be  equal  to  the  number  of  vortex  elements,  and 
would  result  in  the  need  to  solve  large  systems  of  simultaneous  equations. 
For  zero  precision  number,  for  example,  there  would  nave  to  be  800  control 
points  rather  than  80,  and  there  would  be  800  simultaneous  equations  to 
solve. 

We  will,  therefore,  employ  a  modal  approach,  which  greatly  reduces  the 
number  of  unknowns.  The  continuous  distribution  of  circulation  over  the 
span  and  chord  is  assumed  to  be  given  by  a  series  of  known  forms  with 
unknown  coefficients,  following  the  classical  work  of  Glauert  [4].  The  non- 
dimensional  circulation  distribution  due  to  the  angle  of  attack,  and 

/  r\ 

flap  deflection,  Gv  ,  is  assumed  to  be  given  by  the  following  series: 


G(C°  (x,z)  -  T<a’  .  f  Y  c <«>  £„(?)  p  (20  (3.1) 

4U  k-1  &=1 


G«>  (,.,)  .  j[  jx  c<«  £t(J,  ,t(S0  .  (3.2) 

In  the  above  equations  c^are  unknown  mode  amplitudes,  fk(z)  are  bhe 

*\j 

spanwise  modes,  and  P^(s)  are  the.  chordwise  modes.  The  spanwise  modes  are 
given  by  the  following  expression: 


fk(z)  -  sin [(2k  -  l)z]  , 


(3.3) 
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where  z  is  the  transformed  spanwise  coordinate 

r\j  _  * 

z  *  cos  (-2)  . 

The  chordwise  modes  are 


P^)' 


2  r  1  +  cos  s  >. 
- —  ( - ^ — ) 

nc(z)  sin  s 


'u. 


Hfo  -  -4 


TTc(z) 


2*2,  3,  L  -  1 


'V 

PL(t) 


tixt(z) 


'U 

.  1  +  cos  t  . 

'  <\j  ' 

sin  t 


where 


‘V 

t 


*  cos  1[1  - 


*  cos  *[1  - 


2(x  -  xt (z) ) 


2x 


xT(z) 


1  • 


(3.4) 


(3.5) 


(3.6) 


(3.7) 


(3.8) 


(3.9) 


The  last  chordwise  mode,  corresponding  to  2  =  L,  contains  a  square-root 
singularity  at  the  leading  edge  of  the  flap  and  is,  therefore,  omitted  in 
(3.1).  The  first  six  spanwise  and  chordwise  modes  are  plotted  in  Fig.  4.  One 
can  readily  see  from  Fig.  4  that  control  points  must  be  located  very  near  the 
tip  in  order  to  resolve  the  fifth  and  sixth  spanwise  modes. 

The  chordwise  modes  may  be  integrated  to  obtain  the  total  circulation 
around  any  element  of  chord  length: 


^(s)  =  /  p^(s)dx  =  —(s'  +  sin  s) 


_  As  f  As  ,  1A  sin  2ss 

P2(s)  =  /  p2(s)dx  =  ~(s - 2 - 


(3.10) 

(3.11) 
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i  : 


1  H 


1 

'  f 

it 


•i 

* 

■i  \ 


r  _  _  2,sin(£-2)s  sin 

f  p£(s)dx  —(  2(^12) - 2T"> 


£  *  2,  3,  ....  L  -  1 
(3.12) 


PL(t) 


/  pL(t)dx 


1  A/ 

“(t  +  sin  t) 


(3.13) 


The  total  circulation  around  any  spanwise  station  due  to  each  mode  is 
obtained  by  substituting  the  leading  and  trailing  edges  as  limits  in  the 
preceding  integrals.  It,  therefore,  follows  that 


*j(*) 

/  p1(s)dx  -  1 
xL(z) 

^(z) 

\  f  P9(s)dx  =  1 

!  V2> 

I 


xT(z) 

/  p^(s)dx  =  0  £  *  2,  ...»  L  -  1 

*lW 


^(z) 

/  p  (t)dx  =  1  .  (3.14) 

Ju 

0 


We  can  now  combine  (3.14)  with  the  spanwise  modes  to  obtain  an  expression  for 

the  spanwise  circulation  distribution 

xT(z)  K 

Fv^  (z)  =  /  (x,z)dx  =  4U  l  (c£j°  +  ck2^)fk(z).  (3.15) 

Xj^(z)  k=l 


13 


(3.16) 


We  can  now  obtain  the  circulation  strengths  of  each  discrete  spanwise 
element  corresponding  to  each  mode  of  unit  amplitude.  The  circulation  must 
be,  of  course,  constant  over  the  span  of  the  element.  We  chose  to  make  this 
value  correspond  to  that  of  the  continuous  distribution  at  the  midspan  of  the 
element.  Consequently,  the  spanwise  mode  value  is  obtained  by  substituting 

'U 

the  value  of  z  corresponding  to  the  midspan  of  the  element  in  question  into 
(3.3).  The  circulation  is  then  obtained  by  multiplying  this  value  from  (3.5) 
by  the  integral  of  the  particular  chordwise  mode  over  one  chordwise  interval 
from  (3.10)  -  (3.13).  The  sum  of  the  strengths  of  the  spanwise  vortices  over 
any  chordwise  panel  will,  therefore,  be  equal  tc  the  total  circulation  around 
the  midspan  of  the  panel  in  the  continuous  case. 

The  circulation  of  the  twc  trailing  vortices  shed  from  the  ends  of  each 
spanwise  vortex  segnent  must  then  have  the  same  magnitude  with  an  appropriate 
algebraic  sign.  Finally,  since  corresponding  elements  on  the  image  portion 
of  the  lifting  surface  have  the  same  stre^gcn,  this  computation  needs  to  be 
made  only  over  the  semispan. 

An  alternative  procedure  for  determining  individual  vortex  strengths 
is  a  method  originated  by  Faulkner  [1).  The  relative  vortex  strengths 
corresponding  to  any  chordwise  mode  are  obtained  by  requiring  that  the 
downwash  at  each  panel  boundary  be  exact  in  two-dimensional  flow.  The 
results  are  given  in  Table  1  for  8  chordwise  panels,  and  in  Table  2  for 
50  chordwise  panels.  The  only  significant  differences  occur  in  the  first 
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chordwise  mode  near  the  leading  and  trailing  edges,  and  both  methods  result 
in  essentially  identical  final  values  of  lift  distribution.  However ,  the 
Faulkner  method  requires  a*  greater  amount  of  com_  atation,  since  the  relative 
panel  widths  on.  the  skeg  and  flap  are  functions  of  spanwise  position.  For 
zero  precision  number,  the  Faulkner  method  would  require  the  solution  of 
120  different  sets  of  simultaneous  equations,  each  with  50  unknowns.  With¬ 
out  the  complication  of  the  flap.,  the  Faulkner  procedure  for  obtaining 
relative  vortex  strengths  would  be  much  simpler,  since  the  relative  vortex 
strengths  could  be  pre-computed . 
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4.  COMPUTATION  OF  INDUCED  VELOCITIES 


The  lattice  arrangement  as  described  in  Section  2  results  in  a  set  of 
discrete,  skewed  horseshoe  vortex  elements,  as  shown  in  Fig.  2.  Consider 
a  particular  horseshoe  element  consisting  of  a  spanwise  vortex  of  strength 
F  extending  from  (xi,zi)  to  (x2,Z2) ,  joined  by  semi-infinite  trailing  vortices 
of  strength  V  starting  at  (x2,z2)  and  -T  starting  at  (xi,zi).  For  each 
such  element  located  in  the  interval  0  <  z  <  1,  there  will  be  a  corresponding 
image  element  of  the  same  strength  extending  from  (X2,~Z2)  to  (xi,-z'x).  The 
trailing  vortex  shed  at  (x2»-Z2)  will  have  a  strength  ~r,  while  the  vortex 
shed  from  (xx,-zi)  will  have  a  strength  of  +F. 

Let  us  first  consider  the  spanwise  vortex  segment.  Defining  (£,£)  as  the 
coordinates  of  a  general  point  on  the  vortex,  and  (x,z)  as  the  coordinates  of  the 
control  point,  the  velocity  may  be  written  in  accordance  with  the  law  of  Biot- 

(a), 


Savart  as  follows: 

*  r 


(x-g)dC  -  (z~E)d£ 

[(x-V)"2  +  C"z-c>23  3/2  * 


(4.1) 


Xl,Z2 


Along  the  vortex,  we  have: 


t  =  JiL 

dC 


x,-x, 

- *—  ~  const 

z2-zi 


(4.2) 


so  that  (4.1)  may  be  expressed  xn  terms  of  the  variable  £  alone  and  readily 

integrated  to  give  the  result: 

,  \  22 
4 it  v^  (x, z) _ 2a£+b _ 

T  2d/aC2+b£+c  „ 
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where 

a  -  1+t2 
e  *  x~x1+tz1 
b  *  -2(et+z) 
c  *  e  +z 

d  *  e-tz  .  (4.3) 


Equation  (4.3)  is  not  suitable  for  numerical  computation  if  d  becomes 
small.  As  is  evident  from  Fig.  2,  d  is  the  horizontal  distance  between 
the  control  point  and  the  spanwise  vortex  line  or  its  extension.  For  the 
vortex  arrangement  shown  in  Fig.  3,  d  can  never  be  less  than  approximately 
one  half  the  chordwise  grid  spacing  for  vortex  elements  occupying  the  semi¬ 
span  interval  0  <  z  <  1.  The  distance  d  will  only  become  small  enough  to 
cause  problems  if  the  aspect  ratio  is  extremely  high.  For  this  situation, 
the  limiting  form  of  (4.1)  valid  for  d«l  is: 


4ir  v 


(s) 


2 71  fi 
r  Id 


_2sL 


2d 


(2az1+b)2  (2az2+b)' 


(4.4) 


provided  zx  <  z  <  z2.  In  that  case,  the  velocity  approaches  that  of  an 
infinite  vortex  as  given  by  the  first  term  in  (4.4).  The  second  and  third 
terms  represent  small  corrections  to  the  result  for  an  infinite  vortex.  It 
has  been  found  by  numerical  experimentation  that  round-off  error  is  minimized 
if  (4.4)  is  used  for  |d|  <  0.002. 

If  z  is  not  in  the  interval  z\  <  z  <  z2,  the  first  term  of  (4.4) 
disappears,  and  the  velocity  tends  to  zero  as  d  -*■  0.  If  (4.3)  is  used  in 
this  case,  catastrophic  round-off  error  can  occur.  This  is  because  the 
integral  becomes  very  large,  but  Independent  of  ? .  Hence,  the  integral 
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from  zx  to  z2  becomes  the  small  difference  between  the  two  large  numbers. 

This  situation  does  not  occur  in  practice  for  elements  over  the  semi-span. 
However,  due  to  the  inclination  of  the  spanwise  vortex  lines,  it  is  evident 
from  Fig.  3.  that  a  control  point  could  easily  be  aligned  with  the  extension 
of  an  image  vortex  segment  located  in  the  interval  -1  <  z  <  0.  This  can 
result  in  seemingly  random  errors  which  come  and  go  with  minor  changes  in 
grid  spacing  or  planform.  This  problem  is  eliminated  entirely  if  the  velocity 
is  set  to  zero  for  |dj  <  0.002  provided  z  is  outside  the  interval 
zi  <  z  <  Z2.  Inclusion  of  the  small  correction  similar  to  the  second  term 
of  (4.4)  makes  no  difference,  and  is  therefore  an  unnecessary  complication. 

The  velocity  induced  by  the  two  trailing  vortices  is: 


(4.5) 

Since  the  trailing  vortices  are  all  parallel,  the  round-off  error  situation 
is  much  simpler  than  for  the  spanwise  vortices.  Difficulties  could  only  be 
encountered  for  aspect  ratios  far  below  those  of  practical  interest.  Equation 


(4.5)  may  therefore  be  used  for  all  elements. 

The  total  velocity  induced  by  a  horseshoe  element  is  the  sum  of  either 
(A. 3)  or  (A. 4)  and  (4.5).  The  contribution  of  corresponding  image  elements 
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may  then  be  obtained  by  repeating  this  computation  following  the  substitution: 

xi  x2 


z\  -*■  -z2 
X2  -*■  X! 
Z2  -Zi 


(4.6) 


which  results  in  the  correct  algebraic  sign  for  each  of  the  individual  ele¬ 
ments.  We  shall  use  the  notation  v  .  .  to  denote  the  velocity  induced 

m,n,i,j 

at  the  control  point  located  at  (za»xn)  by  the  complete  unit  strength 
horseshoe  element  located  between  (z^x^  j)  and  (z^+^,x_j,+^  ^ )  together 
with  its  image. 


The  FORTRAN  function  HSVEL  listed  in  the  appendix  performs  this  calcu¬ 


lation. 


The  strength  of  the  (i»j)?th  vortex  element  corresponding  to  the  k'th 
spanwise  and  £'th  chordwise  modes  of  unit  amplitude  is 


fk(ziJ  P^(si,j)  » 


(4.7) 


where  z^  is  the  transformed  coordinate  of  the  midspan  of  the  element 


z^  =  cos  (- 


(zi  +  W 


(4.8) 


and  P  ^(si  ^)  is  the  integral  of  the  chord  load  function  obtained  by  substituting 

the  x-coordinates  of  the  leading  and  trailing  edges  of  the  mid span  of  the 

(i,j)'th  element  as  limits  of  integration  in  (3.14). 

The  velocity  induced  at  the  (m,n)'th  control  point  by  the  (k,&)*th  mode  of 

unit  amplitude  is,  therefore,  obtained  by  summing  the  product  of  the  mode 

strengths  from  (4.7)  with  the  velocities  for  unit  circulation  obtained  from 
(4.3)  -  (4.5)  over  all  I  x  j  horseshoe  elements: 


_ 


1  -v  J  ^ 


(4.9) 


Since  the  contribution  of  the  image  is  included  in  v  .  , ,  the  summation  in 

(4.9)  is  only  over  the  elements  in  the  semi-span  interval  0  <  z  <  1. 

We  can  now  write  the  final  expression  for  the  velocity  induced  at  the 
(m,n)'th  control  point  in  terms  of  the  unknown  mode  amplitudes  cjt  ^  : 


K  L 


v  «  y  y  «.  v 

”>n  A  A  k'£ 


(4.10) 


Equation  (4.10) represents  a  set  of  simultaneous  equations  for  the  mode 
amplitude  coefficients,  once  the  values  of  n  are  prescribed  by  the  boundary 
conditions  of  the  problem.  The  mode  amplitude  coefficients  corresponding  to 
zero  flap  deflection  and  unit  angle  of  attack,  c^  are  obtained  by  solving 

(4. 10) with  Vffi  n  =  +1  for  all  values  of  (m,n) .  Similarly,  the  mode  amplitude 
coefficients  for  zero  angle  of  attack  and  unit  flap  deflection  are  obtained 
by  setting  V  =  0  for  all  values  of  (m,:*'  corresponding  to  control  points  on 
the  skeg,  and  V  =  +1  for  all  values  of  (m,n)  corresponding  to  control  points 
on  the  flap. 

If  the  number  of  modes  is  equal  to  the  number  of  control  points,  the 
boundary  condition  may  be  satisfied  exactly  at  all  the  control  points.  If 
the  number  of  modes  is  less  than  this,  (4.10) may  be  solved  by  least  squares 
to  provide  the  closest  possible  fit  at  all  the  control  points.  The  latter 
approach  is  preferred,  since  the  higher  modes  amplitudes  are  generally  of 


.H. A.’/ *►!£?**  . 


questionable  accuracy.  For  this  application  we  solve  (4.10)  for  K»L  unknown 
mode  amplitude  coefficients  by  least  squares  through  80  control  points. 


r 
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5.  COMPUTATION  OF  FORCES 


Lift  and  induced  drag  nay  be  obtained  most  directly  from  the  spanwise 
circulation  distribution  given  in  (3.15)  and  (3.16),  employing  the  well- 
known  results  of  classical  lifting-line  theory  [4] .  The  local  lift  coef¬ 
ficient  per  unit  angle  of  attack  is: 


c  ,Z)  m  =bK(?J . „  _L_  f 

I^i  ^pU2c(z)  c(z)  k£x 


i  *  4a)> 


(5.1) 


noting  that  the  coefficients  c^  and  the  circulation  Y  have  been  obtained 
from  (4.10)  for  an  angle  of  attack  of  unity.  Similarly,  the  lift  coefficient 
due  to  unit  flap  defiction  is: 


CL6(Z>  ■  Sw 


l  Vz>tckf  +  +  ckL>)  •  <5-2) 


Overall  lift  coefficients  are  obtained  by  multiplying  (5.1)  and  (5.2) 
by  the  local  chord,  c(z),  integrating  over  the  span,  and  dividing  by  the  area 
of  the  lifting  surface.  Since  only  the  first  spanwise  mode  contributes,  we 


obtain  the  result: 


-  ,.(■!?  +  ='“>) 
cIS .  +  c<«  +  .<«,  , 


(5.3) 


where  a  is  the  effective  aspect  ratio. 

The  induced  drag  coefficient,  in  accordance  with  lifting  line  theoLy, 
may  be  written  as  follows  in  terms  of  the  present  notation: 
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<L*  K  c£aWo)  2 

iT 11  +  J  {%rV3 

k-2  *U  +C}2' 

ua  U+  z  (2k  1)  .  (5)  (6)  (6)>  j 

cn  12  il 


(5.4, 


The  efficiencies  n  *  of  the  lifting-surface  are  defined  as  the  re¬ 
ciprocals  of  the  quantities  in  square  brackets  in  (5.4),  and  are  equal  to 
one  if  the  spanwise  circulation  coefficients  are  zero  for  k>l. 

To  obtain  the  chordwise  position  of  the  center  of  pressure  at  any  span- 
wise  location,  we  must  obtain  the  first  moments  of  the  chordwise  mode  func¬ 


tions,  p^. 


T 

Mj  »  |  x-p^s)  dx  -  xL(z)  +  ^ 


T 

M2  =  j  x*p2(s)  dx  -  x^Cz)  +  ^2^" 


“.-j*- 


p3(s)  dx  =  - 


Mt  '  \  x- 


p^(s)  dx  =  0 


for  £  =  4 , . . .L-l 


r1  -v,  y*> 

K  “  x.pT(t)  dx  =  — t — 


(5.5) 


a  realistic  approximation  to  the  characterists  of  a  flapped  control  surface 
can  be  made,  provided,  of  course,  that  stall  does  not  occur.  The  constants 
appearing  in  (5,9)  were  obtained  from  experimental  airfoil  data  with  "standard 
roughness"  given  in  [3],  A  sample  tabulation  of  this  type  is  given  in 
Appendix  3, 


6.  TESTS  OF  PROGRAM  ACCURACY 


The  effect  of  the  spanwise  and  chordvise  grid  spacing  on  computed 
lift-slope  for  a  typical  control  surface  is  given  in  Table  3.  Variations 
with  precision  numbers  are  quite  small,  and  it  may  be  concluded  that  zero 
precision  uumbers  are  satisfactory  for  plauform  shapes  similar  to  the  case 
examined . 

Table  3 

Variations  of  Parameters  with  Extreme  Precision  Numbers 


*v 

CDia 

CDi6 

CLa 

CL6 

ua 

6 

0 

0 

0 

0.114 

0.116 

3.138 

1.776 

0.996 

0.984 

0 

2 

0.114 

0.116 

3.132 

1.877 

0.996 

0.980 

2 

0 

0.114 

0.116 

3.101 

1.756 

0.995 

0.980 

a  « 

2.8 

A  * 

15° 

f  ® 

0.2 

X  * 

0.6 

Convergence  of  the  solution  with,  increasing  numbers  of  elements  is  a 
necessary  but  obviously  not  sufficient  test  of  program  accuracy.  Fortunately, 
the  results  for  zero  flap  deflection  may  be  compared  with  corresponding 
solutions  obtained  by  other  current  numerical  lifting-surface  theory  tech¬ 
niques.  A  recent  publication  by  Langan  and  Wang  [5]  is  particularly  helpful 
in  this  regard.  This  reference  compares  the  spanwise  distributions  of  lift 
for  a  tapered  wing  of  aspect  ratio  5  obtained  by  fifteen  different  lifting- 
surface  computer  programs.  Fig.  5  shows  the  results  of  the  present  program 
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plotted  on  a  reproduction  of  Langan  and  Wang's  Fig.  13  vhich 
are  In  excellent  agreement  with  the  average  results  of  the  fifteen  programs 
compared  In  [5].  The  overall  lift  coefficient  obtained  by  the  present 
program  for  this  example  is  4.106  which  agrees  exactly  wic  the  results  given 
in  [5]  obtained  by  the  Tulinius  program!  Induced  drag  was  also  found  to 
be  in  good  agreement. 

A  comparison  was  also  made  of  the  lift  coefficient  of  a  rectangular 
planform  of  unit  aspect  ratio.  Watkins,  Wools ton,  and  Cunningham  in  a 
1959  report  [6],  cite  the  following  values  obtained  from  several  sources, 
to  which  we  have  added  the  results  from  the  present  program: 


Source 

Jones  (low  aspect  ratio  limit) 

1.571 

Lawrence 

1.400 

Hsu 

1.497 

Watkins-Woolston-Cunningham 

1.455 

Present  Program 

1.508 

As  a  further  check  in  the  low  aspect  ratio  range,  J.  Dulmovits  of 
the  Grumman  Aircraft  Engineering  Corporation  kindly  offered  to  run  his 
program  for  a  tapered  planform  of  effective  aspect  ratio  2.8.  He  obtained 


IT 


(6.1) 


0.967 

o 


0.970 

o 


Fig.  5  Spanwise  Lift  Distribution  Obtained  by  Present  Program 
Superimposed  on  Fig.  13  of  [5j 
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SYMBOL 

PROGRAM 

O 

TULINUS 

□ 

DULMOVITS 

0 

MARGASON-  LAMAR 

A 

Gl  ESING 

b 

RUBBERT 

Q 

LOPEZ-SHEN  (EVD) 

0 

H  AVI  LAND 

Q 

JORDAN 

Q 

LAMAR  (MULTHOPP) 

Ci 

WIDNALL 

o 

BANDLER  (ERA) 

Q 

ROWE 

0 

CUNNINGHAM 

0 

JACOBS-TSAKONAS 

0 

LOPEZ  (KUCHE  MANN) 

Table  4  Identification  of  Symbols  in  Fig,  5 
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where  sfi  is  the  transformed  chordwise  coordinate  evaluated  at  the  position 
of  the  flap  hinge.  For  a  50%  flap  the  present  program  gives  C^/C^  *  .807 
compared  with  a  value  of  .818  obtained  from  (6.1).  For  a  20%  flap  the  com¬ 
puted  value  is  .526  vs.  ,550  obtained  from  (6.1).  We  should  expect  slightly 
lower  results  for  an  aspect  ratio  of  60  so  that  this  result  seems  reasonable. 

The  calculation  for  an  aspect  ratio  of  60  also  provided  an  additional 
check  on  in  the  limit  of  high  aspect  ratio.  The  value  obtained  was  5.915 
which  compares  closely  with  a  value  of  6.080  which  one  would  obtain  from 
the  lifting-line  theory. 

It  may  also  be  of  some  interest  to  compare  the  value  of  given 

by  the  present  program  with  that  given  by  Equation  (23a)  of  Chapter  VIII  of 
Principles  of  Naval  Architecture  [8] .  The  latter  is  a  simple  empirical 


equation: 


(0.9)  (2TT)a 

(cosA  /-5LT4)  +1.8 
r  cos*A 


(6.2) 


which,  for  an  aspect  ratio  a  ■  2.8  and  a  sweep  angle  A  «*  15°  gives  a 
value  of  2.997,  which  is  about  five  percent  below  the  computed  value  of 
3.138.  Since  (6.2)  was  developed  to  provide  correlation  with  experimental 
results,  the  lower  value  of  lift  slope  is  not  at  all  unreasonable. 


7.  DESIGN  APPLICATION 


The  present  program  was  initially  developed  to  provide  a  rational 
basis  for  selecting  optimum  planforms  for  an  experimental  series  of  flapped 
rudders  f 92 •  The  first  two  rudders  in  this  series  were  to  have  an  effective 
aspect  ratio  of  2.8  and  flap  area  ratios  of  20%  and  10%,  respectively.  These 
were  planned  as  an  extension  of  an  earlier  series  of  experiments  on  flapped 
rudders  published  in  1972  [10] . 

With  such  small  flaps,  minor  changes  in  sweep  angle  and  taper  ratio 
have  a  large  effect  on  the  spanwise  distribution  of  flap  chord.  The  objective 
is  to  determine  the  sweep  angle  and  taper  ratio  in  such  a  way  as  to  optimise 
the  following  parameters: 

a)  Maximize  lift  slopes  and  ; 

b)  Minimize  induced  drag; 

c)  Provide  nearly  uniform  distribution  of  over  span. 

This  requirement  is  based  on  a  decision  to  use  a  constant 
thickness /chord  ratio  of  15%  over  the  span.  Uniform 
lift  coefficient  is  then  optimum  both  for  delay  of  cav¬ 
itation  inception  and  delay  of  stall. 

d)  Provide  sufficient  flap  tip  chord  to  permit  installation 
of  a  hinge. 

Calculations  were  made  for  five  combinations  of  taper  ratio  and  sweep 
for  the  20%  flap  rudder  and  the  principal  results  are  given  in  Table  5. 

Plots  of  spanwise  distribution  of  lift-slopes  are  shown  in  Fig.  6. 


It  is  clear  from  Table  5  that  overall  lift  and  drag  characteristics 
are  very  insensitive  to  sweep  and  taper  within  the  fairly  limited  range 
permitted  due  to  the  small  flap.  This  is,  of  course,  a  characteristic 
of  low  aspect  ratio  lifting  surfaces,  so  that  this  result  is  not  surprising. 
If  one  looks  closely,  one  can  see  the  trend  of  decreasing  lift  slope  and 
increasing  drag  with  increasing  sweep  angle. 

The  effect  of  taper  ratio  on  spanwise  distribution  of  lift  is  much 
more  pronounced,  as  is  evident  from  Fig.  6.  Here  it  is  clear  that  a  taper 
ratio  of  0.9  unloads  the  tip  too  much,  while  a  taper  ratio  of  0.5  does  the 
reverse,  and  that  a  value  of  0.6  seems  about  right.  A  change  ^.i  sweep  angle 
from  11  to  18  degrees  has  essentially  no  effect  on  CLa  >  and  a  small 
effect  on  . 

Rudder  No.  2  appears  to  have  the  best  spanwise  distribution  of  lift  as 
well  as  maximum  lift  and  minimum  drag.  However,  the  tip  chord  of  the  flap 
is  extremely  small,  which  would  cause  difficulties  in  the  model  and  possibly 
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in  the  full-scale  hinge  design.  An  increase  in  sweep  angle  from  11  to  15 
degrees  overcomes  this  problem  with  very  little  compromise  in  performance. 
Rudder  No.  3  was  therefore  selected  for  the  test  program.  This  has  the 
additional  practical  advantage  of  a  nearly  vertical  trailing  edge  and  a 
constant  flap  chord. 

The  constraint  on  flap  chord  for  the  10%  flap  area  ratio  rudder  is 
even  more  severe.  Rudder  No.  6  shown  in  Table  5  and  Fig.  6  has  the  same 
planform  as  No.  3,  but  with  a  flap  area  of  10%.  Again,  we  find  that  CL$ 
is  a  little  too  low  at  z  »  0.  Increasing  this  by  a  decrease  in  sweep  angle 
is  even  more  impractical  in  this  case  due  to  the  small  flap  chord.  We 
conclude,  therefore,  that  this  planform  seems  to  be  nearly  optimum  for  both 
the  20%  and  10%  flap  rudders. 
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to 
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Fig.  6 


Spanwise  Distribution  of  Lift  due  to  Angle  of  Attack 
and  Flap  Deflection  for  Rudder  Designs  of  Table  6 
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APPENDIX  1 


Instructions  for  Preparing  Computer  Program  Input  Data 


Integer  variables  oust  be  right-adjusted  in  their  fields  -  Real  variables 
must  contain  a  decimal  somewhere  within  their  field. 

A.  FIRST  CARD 


SYMBOL  MODE  FIELD  LIMITATIONS  DESCRIPTION 


RDM 

Integer 

4  |  1-6 

Number  of  spanwise  modes-K 
Reconmended  value  is  6 

LT 

Integer 

8  ;  3 

Number  of  chordwise  modes-L 
Recommended  value  is  6 

IHF 

Integer 

12  0,1  or  2 

! 

Chordwise  precision  index  1^ 
Recommended  value  is  0 

1 

IV 

Integer 

16  0,1  or  2 

1 

Spanwise  precision  index  ! 

Recommended  value  is  0 

ASR 

Real 

17-24  ,  * 

Geometric  aspect  ratio,  |a 

AF 

Real 

25-32  0.1<AF<0.9 

Flap  area  ratio,  f 

T 

Real 

33-40  * 

Taper  ratio,  X  j 

1 

PP 

Real 

41-48  * 

Sweep  angle  in  degrees,  A 

KOPT 

Integer 

56  0  or  1 

i 

1 

i 

i 

# 

* 

\ . . .  . . 

If  K0PT=1,  the  subroutine  • 

OPTION  will  be  called  and 
i  will  perform  the  calcula¬ 
tions  appearing  on  the  third  1 

page  of  the  output 
_ 

B.  SECOND  CARD 

This  card  contains  the  indices  of  the  ten  chordwise  panels  which  contain 
control  points,  NCP(N) ,N=1,10.  These  are  integers  which  are  read  in  ten 


Combinations  of  these  quantities  must  be  so  chosen  that  flap  chords  are 
positive  for  all  z. 
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fields  of  8  columns.  The  recommended  values  for  zero  precision  number 
are  3,  8,  13,  18,  23,  28,  33,  38,  42,  and  48.  However,  the  intention  is 
that  these  may  be  altered  for  a  more  advantageous  placement  relative  to 
the  flap  hinge  if  so  desired. 

C.  MULTIPLE  RUNS 

Upon  completion  of  the  calculation,  the  program  returns  to  the  read 
statement  for  tn>  l'rst  card.  A  blank  card  (or  more  specifically,  a  zero 
value  for  KDM)  terminates  the  run. 
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Timing  Information:  0  precision  with  36  modes; 

1.081  minutes  execution  time. 

IBM  System  370/165 

Memory  Requirement  100K 
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APPENDIX  4 


Control  points  should  be  located  at  a  point  within  the  element  where  the 
effects  of  the  velocities  induced  by  all  four  boundaries  of  the  element  cancel. 
This  is  essential  in  order  for  the  discrete  system  to  converge  to  the  Cauchy- 
Principal  Value  of  the  corresponding  continuous  singular  integral.  For  a  rec¬ 
tangular  element  this  will  be,  of  course ;  at  the  centroid  of  the  element. 

With  the  present  vortex  scheme  the  elements  may  be  both  swept  and  tapered, 
so  that  a  question  arises  as  to  whether  the  placement  of  the  control  point  is 
critical.  A  detailed  calculation  was  therefore  made  for  an  extreme  element 
with  geometric  characteristics  as  shown  in  Fig.  7. 

The  velocity  induced  by  the  four  line  vortex  elements  was  computed  at 
finely  spaced  intervals  throughout  the  interior  of  the  element.  In  this  case, 
the  null  point,  0,  was  found  to  be  displaced  from  the  mid-chord /mid-span  point, 
M,  by  the  amount  shown  in  Fig.  7.  This  deviation,  as  well  as  the  velocity  at 
M  is  extremely  small. 

In  order  to  obtain  an  estimate  of  the  error  introduced  by  locating  control 
points  at  M  rather  than  at  0  in  the  complete  lifting  surface  program  a  test 
run  was  made  with  all  control  points  displaced  by  two  percent  of  their  panel 
chords.  This  resulted  in  a  change  of  0.92%  in  predicted  lift  slope.  Since 
this  displacement  was  far  in  excess  of  the  value  shown  in  Fig.  7,  it  can  be 
concluded  that  the  error  introduced  by  locating  control  points  at  the  mid-chord 
/mid-span  position  of  an  element  is  negligible. 


Fig.  7  The  Element  Used  for  Estimation  of  the  Error  Caused 
by  Control  Point  Position  Deviations 
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